Nonlocal Form of the Rapid Pressure-Strain Correlation in Turbulent Flows 
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A new fundamentally-based formulation of nonlocal effects in the rapid pressure-strain correlation 
in turbulent flows has been obtained. The resulting explicit form for the rapid pressure-strain 
correlation accounts for nonlocal effects produced by spatial variations in the mean-flow velocity 
gradients, and is derived through Taylor expansion of the mean velocity gradients appearing in 
the exact integral relation for the rapid pressure-strain correlation. The integrals in the resulting 
series expansion are solved for high- and low-Reynolds number forms of the longitudinal correlation 
function f(r), and the resulting nonlocal rapid pressure-strain correlation is expressed as an infinite 
series in terms of Laplacians of the mean strain rate tensor. The new formulation is used to obtain 
a nonlocal transport equation for the turbulence anisotropy that is expected to provide improved 
predictions of the anisotropy in strongly inhomogeneous flows. 
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the pressure-strain correlation tensor 
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By far the most practical approaches for simulat- 
ing turbulent flows are based on the ensemble-averaged 
Navier-Stokes equations. However, such approaches re- 
quire a suitably accurate closure model for the Reynolds 
stress anisotropy tensor ay, defined as 



n y (x) = -p'(x)^,(x). 



(i) 
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where u^u'a are the Reynolds stresses and k 
the turbulence kinetic energy. Over the past half cen- 
tury, a wide range of closures for ay have been proposed. 
Of these, so-called Reynolds stress transport models that 
solve the full set of coupled partial differential equations 
for ay are currently regarded as having the highest fi- 
delity among practical closures. Such closures start from 
the exact transport equation for ay, namely 
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where for clarity we are restricting the presentation to 
incompressible flows. In (2), D/Dt is the mean- flow ma- 
terial derivative, fy = — (w-uj duj/dxi + u'jU[ dui/dxi) 
is the production tensor, ey is the dissipation tensor, and 
all remaining viscous, turbulent, and pressure transport 
terms are contained in Dy, with P = P nn /2, e = e nn /2, 
and D = D nn /2. In such Reynolds stress transport clo- 



sures, the production tensor needs no modeling since u^Uj 
is obtained from ay , and standard models for ey and Dy 
are discussed in Rcfs. [1, 2]. The principal remaining dif- 
ficulty is in accurately representing Ily in (2), namely 



where 
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are the strain rate fluctuations. The pressure-strain cor- 
relation has received considerable attention, however de- 
veloping a fundamentally-based yet practically implc- 
mentable form for Ily remains one of the primary chal- 
lenges in turbulence research. 

The difficulty in representing fly stems in large part 
from the inherently nonlocal nature of the pressure-strain 
correlation, since the local pressure p' in (3) depends on 
an integral over the entire spatial domain of the flow. 
Some progress has been made by splitting fly into the 
sum of "slow" and "rapid" parts [3], where the rapid 

(r) 

part Ily is so named due to its direct dependence on 
the mean- flow velocity gradients dui/dxj, variations in 
which have an immediate effect on Daij/Dt. Typically, 

the slow part ETV.-' is represented in terms of the local 
values of ay and e. For the rapid part, it has been com- 
mon (e.g., [3-5]) to take the mean velocity gradients as 
being sufficiently homogeneous that they can be brought 
outside the integral. Under certain conditions [5] the re- 
maining integral can then be solved for the local part of 

(r) 

Ily . This is then typically combined with additional ad 
hoc terms involving ay to model the rapid part solely in 
terms of local flow variables. Together with the assumed 
local representation for the slow part, this yields a purely 
local formulation for Ily that allows (2) to be solved, but 
that neglects all nonlocal effects in the evolution of the 
anisotropy. 

Such purely local models for Ily have allowed rel- 
atively accurate simulations of homogeneous turbulent 
flows, where by construction there are no spatial varia- 
tions in dui/dxj and thereby all nonlocal effects vanish. 
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However most practical situations involve strongly inho- 
mogeneous flows, where large-scale structure and other 
manifestations of spatial variations in the mean-flow ve- 
locity gradients can produce significant nonlocal effects 
in the turbulence, the neglect of which in 11^ can lead 
to substantial inaccuracies in the resulting anisotropy. 
Such nonlocal effects arc significant even in free shear 
flows such as jets, wakes, and mixing layers, and can be- 
come especially important in near- wall flows, where flow 
properties vary rapidly in the wall-normal direction. Im- 
proving the fidelity of turbulent flow simulations requires 
a fundamentally-based formulation for nonlocal effects in 

(r) 

II \- to account for spatial variations of velocity gradients 
in the ensemble-averaged flow. 



Various methods for addressing such spatial variations 
have been proposed, however nearly all suffer from a 
lack of systematic physical and mathematical justifica- 
tion. For near-wall flows, by far the most common yet 
also least satisfying approach is the use of empirical "wall 
damping functions" {e.g., [2]). Although such functions 
are relatively straightforward to implement, they are also 
distinctly ad hoc and as a consequence do not perform 
well across a wide range of flows. Moreover, wall func- 
tions typically conflate the treatment of a number of near- 
wall effects that in fact originate from distinctly different 
physical mechanisms, including low Reynolds number ef- 
fects, large strain effects, and wall-induced kinematic ef- 
fects, and are not formulated to specifically account for 
nonlocality due to spatial variations in the mean flow 
gradients. 



In the following we depart from these prior approaches 
by systematically deriving a new nonlocal formulation 
for the rapid pressure-strain correlation from the exact 
integral relation for the rapid part of Hy. Specifically, 
nonlocal effects due to mean-flow velocity gradients are 
accounted for through Taylor expansion of dut/dxi in 
the rapid pressure-strain integral. The resulting nonlocal 
form of the rapid pressure-strain correlation 11^ appears 
as a series of Laplacians of the mean strain rate tensor. 
The only approximation involved - beyond the central 
hypothesis on which the present formulation is based - 
is an explicit form for the longitudinal correlation func- 
tion /(r), though the effect of this is only to determine 
specific values of the coefficients in an otherwise funda- 

(r) 

mental result for the nonlocal effects in II; - ' . The coeffi- 
cients are obtained here for the exponential form of f(r) 
appropriate for high Reynolds numbers, and for the ex- 
act Gaussian /(r) that applies at low Reynolds numbers. 
The resulting formulation for the rapid part of 11^ then 
provides a new nonlocal anisotropy transport equation 
that can be used with any number of closure approaches 
for representing a,ij , including Reynolds stress transport 
models as well as explicit stress models suitable for two- 
equation closures. 



II. NONLOCAL FORMULATION FOR THE 
PRESSURE-STRAIN CORRELATION 

The starting point for developing a fundamentally- 
based representation for Iljj is the exact Poisson equation 
for the pressure fluctuations p' appearing in (3), namely 



1 9 / „du k du', d 2 
-VV = -2- ' 



dxi dxk dxkdxi 



(uM - u' k u[) (5) 



(e.g., [6]). Beginning with Chou [3], it has been common 
to write p' in terms of rapid, slow, and wall parts as 



p' = p'W + p'^ + p'M 



(6) 



defined by their respective Poisson equations from (5) as 

(7) 



I v y (r) = _2^M 

p dxi dx k 



-W (s) - 
p 



d 2 



dxkdx t 
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u k u i 



i v y (tu) = 0. 
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The effect of p' {w ^ is si gnificant in (2) only in the ex- 
treme near- wall region of wall-bounded flows [6, 7]. The 
remaining rapid and slow parts produce corresponding 
rapid and slow contributions to the pressure-strain cor- 
relation Hij in (3), with Green's function solutions of (7) 
and (8) giving these as 



nS } (*) = z 



du k (x) <9uJ(x) d 3 it 



7r J R dxi dx k 



9 2 K<k a, 

dx k dxt ijK >\*-x 



(10) 



(11) 



where the integration spans the entire flow domain R. 

(s) 

The slow part II \-' is typically not treated in a sys- 
tematic fashion via integration of (11). Instead, nearly 

(s) 

all existing representations for II ^ ■ are based on insights 
obtained from the return to isotropy of various forms 
of initially-strained grid turbulence. The most common 
representation for IT^ is Rotta's [4] linear "return-to- 
isotropy" form 



l4 s) - -deaij , 



(12) 



where all variables are local and C\ is typically in the 
range 1.5 — 1.8 (e.g., [2, 6]). Sarkar and Speziale [8, 9] 
have argued that additional quadratic terms should be 
included in (12), but it has been noted [2] that these 

(s) 

are typically small. As a result, representations for IQ. 
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remain relatively simple, and the form in (12) continues 
to be widely used. 

(r) 

By contrast, IL; ■ has received substantially greater at- 
tention. The direct effect of the mean velocity gradients 
dUk I dxi on this rapid part of the pressure-strain corre- 
lation is apparent in (10). In the following sections, we 
use the integral in (10) to develop a fundamentally-based 

(r) 

representation for II > . ' that accounts for nonlocal effects 
resulting from spatial nonuniformities in the mean veloc- 
ity gradients. 



where k is the turbulence kinetic energy, it can be shown 
[5] that Mujk in (16) becomes 



Mujk = Y^fc (45j k Sii - 5 i3 5 k i - SjiSki) , (20) 

where the leading k again denotes the turbulence kinetic 
energy. Using (20) in (17) then gives the rapid pressure- 
strain correlation as 



(21) 



A. Prior Local Formulation for n^(x) 

Chou [3] first suggested the notion of using the inte- 
gral form in (10) to obtain a representation for the rapid 
pressure-strain correlation. Subsequently, Rotta [4] and 
then Crow [5] used that approach to rigorously derive the 



where Sij is the local mean-flow strain rate tensor 



- If d%h du 3 

on 



2 \dxj dxi 



(22) 



Typically, (21) is used as the leading-order isotropic 
term in tensorial expansions for the rapid pressure-strain 
purely local part of U^> , by assuming the mean veloc- corre l a tion Iljftx), where the remaining terms are ex- 



ity gradients in (10) to vary sufficiently slowly that they 
could be taken as constant over the length scal e on which 
the two-point correlation [du' k (k)/dxi] S'^ (x) in (10) is 
nonzero. Under such conditions, the mean velocity gra- 
dient in (10) can be taken outside the integral, and IT^ 
then becomes 



dxi 



R 



dx k 



With S' i j(x) in (3), the integrand in (13) involves two- 
point correlations among velocity gradients of the form 



9^(x)^;(x) = d 2 R u (r) 
dxj dik drjdrk 



where Ru(r) denotes the velocity fluctuation correlation 



J2«(r) = «{(x)«{(r + x) 
with r = x x. Defining [3-5] 

1 f d 2 R ll (r)d 3 r 



M, 



iljk 



2ir ./r drjdru r 



the rapid pressure-strain correlation in (13) can then be 
expressed as 



n%\x)n^^[M iljk +M jlik ] 



<9ufc(x) 



pressed in terms of the local anisotropy and the local 
mean velocity gradient tensor. Note however that such 
representations are still purely local, since in going from 
(10) to (13) all spatial variations in the mean velocity gra- 
dients duk/dxi over the length scale on which the two- 
point correlations [du' k (x) / dxi] S^ -(x) are nonzero were 
ignored. The resulting neglect of nonlocal contributions 

(13) to from that approximation can lead to substantial 
inaccuracies in many turbulent flows, including free shear 
flows and wall-bounded flows. For example, Bradshaw et 
al. [10] showed using DNS of fully-developed turbulent 
channel flow [11] that the homogeneity approximation 
used to obtain (17) is invalid for y + < 30. It can be 
further shown (e.g., [11, 12]) that the dominant compo- 
nent S12 of the mean strain begins to vary dramatically 
at locations as far from the wall as y + w 60. Compara- 
ble variations in mean velocity gradients are also found 
in turbulent jets, wakes, and mixing layers, where there 
are substantial spatial variations in S , 12 across the flow. 
Indeed in most turbulent flows of practical interest, there 
are significant variations in the mean flow velocity gradi- 
ents that will produce nonlocal contributions to the rapid 
pressure-strain correlation via (10). In such situations, it 
may be essential to account for these nonlocal effects in 

(r) 

H\j to obtain accurate results from any closures based 
on (2). 



(14) 



(15) 



(16) 



(17) 



Using the homogeneous isotropic form of Ru(r), 
namely 



Ru(r) = 3^ 



with 



fir) 



3 w'(x + r)u'(x) 



k 



(18) 



(19) 



B. Present Nonlocal Formulation for n^'(x) 
In the following, nonlocal effects due to spatial varia- 

(r) 

tions in the mean flow are accounted for in IT^ through 
Taylor expansion of the mean velocity gradients appear- 
ing in (10). The central hypothesis in the approach devel- 

(r) 

oped here is that the nonlocality in IIV . is substantially 
due to spatial variations in du k /dxi in (10), and that in 
order to address this effect all other factors in (10) can be 
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adequately represented by their homogeneous isotropic 
forms. This allows a formulation of the rapid pressure- 
strain correlation analogous to that in (21), but goes be- 
yond a purely local formulation to take into account the 
effects of spatial variations in the mean flow gradients. 

We begin by defining the ensemble-averaged velocity 
gradients 



A M = duk/dxi 



(23) 



and account for spatial variations in A k i(x) in (10) via 
its local Taylor expansion about the point x as 



A /M A I \ , dA ki , r mr p d 2 A kl 
A kl (x) = A kl (x) + r m + 



2 dXfyj dxr. 



(24) 



+ 1 (r m r p ) 

7T/\ Qx it-,-, QXiri • 



where r = x x and all derivatives of A k i are evaluated 
at x, and where n is the order of the expansion. As 
n —> oo, the expansion provides an exact representation 
of all spatial variations in A k i(r + x) from purely local 
information at x. Substituting (24) into (10) then gives 



nl; } (x) = E 



where 



(mp...) lvl Ujk 



n=0 



OXrndX'n . . ■ 



(m P ...)Mg + (mp ...) M% 
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'm'p • • • 
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(25) 



<9r 7 r fe 



(26) 

The nth-order term in (25) involves n derivatives of A k i 
as well as n total indices (mp . . .) in 



{mp...) lvl Ujk- 

From the central hypothesis on which the present treat- 

(r) 

ment of nonlocal effects in rQ - is based, we represent 
Ru(r) in (26) by the form in (18). With the relations 



dr 
drj 



JL 
r 



drj _ 



(27) 



' 3 ' Wl 3 

the double derivative of Ru(r) in (26) is then given by 

d 2 R a {v) = k 
3 



drjdr k 



ldf d 2 f d 3 f 



where we have introduced the compact notation 



(28) 



Cijki = SSjkSii - S i:j S k i - SjiS ki - 3a jk 6u + S lj ai k (29a) 
+Sjia lk + S lk a>ij + Shiaij + S kj au - 3(3ujk , 

bijki = SaSjk + 3&jk5il - S i:j ai k - 5i 3 a lk - S^aij (29b) 
-dkiaij - Skjau + Z(3u 3 k , 



Cijkl — Sudj k Piljk ■ 



with 



_ ' *' 3 a 
O-ij = —p^- , Pijkl 



nrjr k ri 



(29c) 
(30) 



Writing the differential in (26) in spherical coordinates as 
d 3 r = r 2 dr dfl, where dfl = sm9 d9 d(j> and r = [0, co), 
6 = [0, 7r], and <fi — [0, 2tt), since /(r) has no dependence 
on 9 or 4> and since dijki, bijki, and Cijki in (29) have no 
dependence on r, the integrals over these terms in (26) 
can be considered separately. Using (28), the integral in 
(26) can then be written as 



(mp.. 



)M, 



(n) 
iljk 



67m! 

+ 



f0 ° ndf , 

r —dr 



+ 



f 

Jo 



Jo ' dr 



dO^i) 



- dn 



dr 3 



where k in the leading factor is the turbulence kinetic en- 
ergy. With the corresponding expression for (mp...)Mj^ k , 
(25) and (31) provide a nonlocal form for the rapid 

(r) 

pressure-strain rate correlation IT^- in terms of the lon- 
gitudinal correlation f(r). 



C. Representing the Longitudinal Correlation f(r) 

As will be seen later, in (31) the integrals over dfi can 
be readily evaluated. Moreover, for n — the integrals 
over dr are independent of f(r), and thus M u °- k can be 
obtained from the general properties 



A : 



Jo 



f(r)dr, /(0) = 1, /(co) = 0. (32) 



However for n > 0, evaluating the integrals over dr to 
obtain ( mp ...)M^ k requires an explicit form for the lon- 
gitudinal correlation function f(r). We can anticipate, 
however, that the precise form may not be of central im- 
portance to our eventual result for since the only 
role of f(r) is to weight the contributions from velocity 
gradients A k i (x + r) around the local point x. It is thus 
likely that the integral scale A in (32) plays the most 
essential role, since it determines the size of the region 
around x from which nonlocal contributions to the in- 

(r) 

tegral for IIV - will be significant. When r is scaled by 
A, the precise form of /(r/A) is likely to be far less im- 
portant for most reasonable forms that satisfy the con- 
straints in (32). 

Despite its fundamental significance in turbulence the- 
ory, the form of f(r) for any r and all Reynolds num- 
bers Re\ has yet to be determined even for homogeneous 
isotropic turbulence. Perhaps the most widely-accepted 
representation for f(r) comes from Kolmogorov's 1941 
universal equilibrium hypotheses. For large values of 
Re/v = k x l 2 kjv and inertial range separations \ v <C r <C 
A, where \ v ~ (^ 3 /e) 1 / 4 is the viscous diffusion scale and 
A is the integral length scale in (32), the mean-square 
velocity difference is taken to depend solely on r and 



5 



the turbulent dissipation rate e, and thus on dimensional 
grounds must scale as 



[u'(x + r)-u'(x)] ^e^V^. (33) 
Expanding the left-hand side of (33) and using (19) gives 
4, 



;k[l-f(r)} 



e 2/3 r 2/3 



(34) 



Defining the proportionality constant in (34) as Cf and 
rearranging gives the inertial range form of f(r) as 



k 3/2/ £ 



2/3 



(35) 



From Hinze [13], a value for Cf can be obtained in terms 
of the Kolmogorov constant K = (8/9a) 2 / 3 f» 1.7, where 
a ps 0.405, as 



81 

C f = _r(4/3)Jf« 2.24 , 



(36) 



where we have used T(4/3) w 0.893. Expressing A in 
terms of fc and e on dimensional grounds as 



A = C A 



fc 3/2 



(37) 



where Ca is a presumably universal constant, then allows 
the inertial-range form of f(r) in (35) to be given as 



/(-/A) = 1 - \CfCT Q 



\ 2 /3 



(38) 



However the form for f(r) in (38) is valid only for 
inertial-range r values, namely A,; « r « A and thus 
for Re A 3 ^ 4 <C (r/A) *C 1. As a consequence, this form 
cannot be used directly to evaluate the r-integrals in (31). 
However, experimental data from a wide range of turbu- 
lent free shear flows (e.g., [14, 15]) and direct numerical 
simulation results for wall-bounded turbulent flows (e.g., 
[11, 12]) show that f(r) can be reasonably represented 
by the exponential form 



/(r/A) = e 



(39) 



as can be seen in Figures l(a)-(c). Moreover, C\ in 
(37) can be chosen to closely match f(r) in (39) with 
the fundamentally- rooted inertial-range form in (38). In- 
deed, Figure 2 shows that with 



C A « 0.23 



(40) 



the exponential form in (39) gives reasonable agreement 
with the inertial-range form in (38) up to r/A 1. This 
exponential form is thus here taken to represent f(r) in 
high-i?eA turbulent flows, and will be used in (31) to 
obtain an explicit form for the nonlocal rapid pressure- 
strain correlation. Since (25) with (31) is a rigorous for- 

(r) 

mulation for within the central hypothesis on which 
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FIG. 1: Comparison of exponential /(r) in (39) with exper- 
imental data from axisymmetric turbulent jet [14] (a) and 
planar turbulent mixing layer [15] (b), and with DNS data 
from turbulent channel flow at Re T — 650 [12] (c). 



the present approach is based, the exponential represen- 
tation for f(r) is the principal additional approximation 
that will be used below in deriving the present result for 
the rapid pressure-strain correlation. 

While the exponential /(r) appears appropriate for 
high Re a, in the Re a — * limit the Karman-Haworth 
equation [16] allows a solution for f(r). Batchelor and 
Townsend [17] showed that when inertial effects can be 
neglected, this equation can be solved exactly, giving a 
Gaussian form for f(r) as 



/(r/A) = exp 



4 / r 
n VA 



(41) 
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FIG. 2: Comparison of inertial-range and exponential forms 
for /(r/A) in (38) and (39), respectively. Note that C\ = 0.23 
in (40) gives reasonable agreement between the two forms in 
the inertial range Re- 3/i < (r/A) « 1 



Ristorcelli [18] has proposed a blended form for f(r) 
that satisfies various conditions placed on /(r), includ- 
ing those in (32), while recovering the Gaussian f(r) in 
(41) as ReA — * and the exponential f(r) in (39) as 
Re \ — + oo. It should be possible to use such blended 
forms for f(r) to obtain a nonlocal pressure-strain cor- 
relation valid for all Reynolds numbers, following the 
procedure developed herein. In the following we obtain 
the nonlocal pressure-strain correlation using the high- 
Reynolds number exponential form in (39), which should 
be accurate for the vast majority of turbulent flow prob- 
lems, and then show how this result can be extended to 
the low- Reynolds number limit using (41). 



D. Resulting Nonlocal Pressure-Strain Correlation 

Using (39), it can be shown that the integrals over dr 
in (31) give 



r n ^-dr = -n\A n , 
dr 



poo j2 f 

l r n+ ^dr=(n + l)\A\ 



[ r n+2 %tdr = -(n + 2)\A n 
Jo dr 

With these results, (31) is then written as 



M ( n ) 
(mp...) lvl Ujk 



67T 
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TmTp • • ■ 
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(42a) 
(42b) 
(42c) 

(43) 



[dijki - (n + l)hjki + (n + 2)(n + V/Cijki] dft ■ 



The remaining integrals over dQ are all of the form 
( r m r p ■ ■ O/ 7 *™ an d can be solved using the general integral 
relations 



/ 



■» y» y» y» 

m 1 p> q> s 



-dCl = 



-dn = o, 



odd 



4tt 



(n + 1) 



M 



(44a) 



(44b) 



' \&mp3qs 



~t~ ^mq^rs • • • ~\~ ' ' ' 

where the double factorial is defined as 
nil = n(n — 2)(n - 4) ■ • ■ 



(45) 



with 0!! = 1 and (—1)!! = 1, and the terms in brackets on 
the right-hand side of (44b) represent all possible combi- 
nations of delta functions for the indices (m,p, q,s,...). 
For any n, there are (n — 1)!! such delta function terms, 
and each term consists of (n/2) delta functions. 

In (43), for n = it can be shown using (44b) that 



M 



(0) • u 

iji 18 s iven b y 

(o) , 2 



M 



rC 1 ) 



Sijdki - SjiSki) 



(46) 



For n = 1, from (44a) mM^' k — 0, as applies to all odd-n 



- (2) 

cases. For n = 2, from (44b) mpM^L is given by 



mp lyJ -Ujk 



M {2) 



(47) 



2A 2 

k-^r[4:5jk5ii6 m p — 3 (SijSkiSmp + SjiSikSmp) 

—24 (SuSjmSkp + SuSkmSjp) + 4:(SijSi m Sk p + Sij5k m Si p 
+Sji5i m Sk P + SjiSkmSip + SikSi m Sjp + 5ikSj m Sip 
+SkiSi m 6j P + SkiSj m Si P + SjkSi, n Spi + SjkSimdip)} . 

Contracting (46) and (47) with A^i and its derivatives as 
in (25) then gives 



M, 



(0) 



Ujk 



M 



(0) 



and 



d 2 A k 



kl 



(2) 



™P ly± iljk ~<~rn,p lvl jlik 



M 



jlik 



(2) 



(48) 



68 
315 



feA V 5ij , (49) 



where we have used Akk = 0. From (48) and (49), 
the first two terms in the present formulation for the 
pressure-strain correlation in (25) are thus given by 



^ n W(x) = ^ + ||A 2 V 2 ^ + . 



1' 



(50) 



The first term on the right in (50) is the same as that 
in (21) obtained by Crow [5] assuming spatially uniform 
mean velocity gradients. Thus the second term in (50) is 
the first-order nonlocal correction accounting for spatial 
variations in the mean velocity gradient field. 
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To obtain the remaining higher-order nonlocal correc- 
tions in (50), it is helpful to contract (43) with the deriva- 
tives of Aki and again use A kk = 0. It is then readily 
shown that all terms involving Ski, <5fcm, Sk P , ■■■ from 
the integral over dfi are zero when contracted with the 
derivatives of Am, and as a result the coefficients in (29) 
can be simplified as 



a-ijki = 45jk6u - SjiSki - b ijk i 



hjki — SuSjk + b* jkl , 



(51a) 



(51b) 



b ijki = 3 a jk$u - SijOiik - Sikaij - 5k 3 au + 3(3u jk , (51c) 



Cijki = ducxjk - Pujk , 



(51d) 



where b* ]kl has been introduced to simplify the notation. 
Using (51) and contracting (43) with the derivatives of 
Aki we thus obtain 



d n Ai 



OXyfiOXp . . . 



{mp...) lvl Ujk 
1 



(52) 



kA n - 



\ d n A kl 


J 


TmTpfq • • • 


OXyfiOXp . . . 


Jn 





Qtt 

•[(2 - n) 5 a 8 3k - 5 ik 8 3 i - (n + 2) b* jkl 
+(n + 2)(n+l)c ijk i}dn. 

From (44a) all odd-n terms in (52) are zero. For even-n, 
the integrals over dfl are readily evaluated using (44b), 
and it can then be shown that (52) becomes 



d n A 



kl 



0x m 0xp . . . 



(mp...) IVI Ujk 



(53) 



(n + 2)(n + 4)(n + 



(n + 5)!! 



^(v 2 )" /2 [^, + ^,]}. 



Adding the corresponding result for M^ k to (53) then 
gives 11^ from (25) as 



1. 



(r) 



oo 



n—0,even 

where the coefficients are 

(„) _ 4(n 2 + 2n + 9) 



a 



3(n + 5)(n + 3)(n+ 1) 



(54) 



(55) 



Since the indices in (55) and (54) are required to be 
even, we can change the index n to (2n — 2) , where then 



n = 1,2,3,.... This gives the final result for the non- 
local rapid pressure-strain correlation from the present 
approach as 



+ 



OO 

n=2 



C (n) A 2n-2 



(v 2 )' 



with 



(n) _ 16n 2 -16n + 36 



3(2n + 3)(4n 2 - 1) 



(56) 



(57) 



where A in (56) is from (37) and (40). In (57) it may 

be readily verified that cf } = 4/5 and cf } = 68/315, 
consistent with (50) and (21). The first term on the right 

in (57) accounts for purely local effects on II , while the 
series term accounts for nonlocal effects. 

The result in (56) and (57) is the first rigorous formu- 

(r) 

lation for the rapid pressure-strain correlation II \-' that 
accounts for nonlocal effects due to spatial variations in 
the mean velocity gradients. Within the central hypoth- 
esis on which the present approach is based, the principal 
approximation used in deriving (56) and (57) is the expo- 
nential form of f(r) in (39) for high-i?eA turbulent flows. 
However, the only effect of this choice of f(r) is in the re- 
sulting coefficients in (57). All other aspects of (56) 
are unaffected by the particular form of f(r), and instead 
result directly from the fundamental approach taken here 
in solving (10) via Taylor expansion of the mean velocity 

gradients duk /dxi to account for nonlocal effects in . 

The coefficients in (57) from the exponential rep- 
resentation of f(r) are shown in Figure 3. It is appar- 
ent that the n = 1 term in (56), which accounts for the 
purely local contribution to II^ as verified in (50), is by 
far the dominant coefficient. The remaining coefficients 
for n = 2,3,4,... correspond to the nonlocal contribu- 

(r) 

tions to H\j ' , and can be seen in Fig. 3 to decrease only 

slowly with increasing order n. However, while C'^ is 
clearly the dominant coefficient, in (56) the remaining 
coefficients are multiplied by successively higher-order 
Laplacians of the mean strain rate field, and thus may 

(r) 

produce net contributions to IIV - that are comparable 
to, or possibly even larger than, the local term due to 
n = 1. 



E. Corresponding Coefficients for ReA — > 

While the coefficients in (57) arc appropriate for 
Re \ 3> 1, in this Section we use the exact Gaussian form 
for /(r) in (41) that applies in the ReA — * limit to 

(r) 

obtain the result for HQ - applicable to low-i?eA flows, as 
may occur in the near- wall region of wall-bounded turbu- 
lent flows. Using this form for f(r), it can be shown that 
for even-n, which are the only nonzero terms from (43) 
due to (44a), the r-integrals in (31) are modified only 
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-H — Exponential /(r) (High Re^) 
-Q — Gaussian f(r) (Low ife^) 




(n) 

FIG. 3: Comparison of rapid pressure-strain coefficients Cv, 
from (57) for the ReA 1 exponential f(r) in (39), and from 
(58) for the ReA — > Gaussian f(r) in (41). 



to obtain certain components of mpMujk in the high- 
and low- Reynolds number limits, respectively. Note that 
Rotta expressed his results [4] in terms of the transverse 
integral scale L rather than the longitudinal integral scale 
A. The two length scales are related by L = 0.5A for in- 
compressible flows, and this relation can be used directly 
in the ReA — ► limit to compare the present result for 

mp^ujk w ^ tne components obtained by Rotta. For 
Re \ 3> 1, the differences between the inertial-range form 
of f(r) used by Rotta and the exponential form in (39) 
used herein give his L as L — 0.43A. Using these rela- 
tions it may be verified that the limited components of 

(2) 

mp^iijk given by Rotta are in agreement with the com- 
plete result in (47) for Re a ^> 1, and with the result for 
i?eA — > when the factor of 2/w is accounted for as noted 
in Section II E. 

The agreement with those components of 
reported by Rotta [4] provides partial validation of 
the present results. However, the present results go 
much further by addressing the complete components 



M (2) 

mp lvl iljk 



of 



(mp...)M iljk 



for all ft, thereby allowing the first com- 



plete formulation of nonlocal effects in the rapid pressure- 

(r) 

by multiplying the previous results in (42) by the factor strain correlation LT^- due to spatial variations in the 
[(2/Vtt)™ (n/2)!/n!] . The result for M$ k in (46) is inde- mean-flow gradients dujdxj. 



iljk 

pendent of the form of f(r) and thus is unchanged in this 

(2) 

limit, but now m P M- lik in (47) is reduced by the factor 



(n) 
iljk ' 



2/tt. With the remaining higher-order terms ( mp ...)M, 

it may be readily verified that the result for Tl\^ in (56) 
is unchanged in this low-i?eA limit, but the coefficients 
C2 are now given by 



16n 2 - 16n + 36 
3(2n + 3)(4n 2 - 1) 



(2n-2)! 



n-l 



(58) 



III. NONLOCAL ANISOTROPY TRANSPORT 
EQUATION 

The present result for nonlocal effects in the rapid part 
n\j of the pressure-strain correlation, given by (56) with 

the coefficients C ( 2 n) in (57) or (58) and with A in (37) 
and (40), can be combined with (12) for the slow part 
in (2) as 



n< s) to give n y 



where again n = 1,2,3,.... The effect of the additional 
factor in (58) relative to (57) is to damp the higher-order 
terms in the ReA 3> 1 coefficients, as shown in Fig. 3. It 
is apparent that in this ReA —> limit, only the first non- 
local term (n — 2) in (56) is significant, with all higher- 
order coefficients being negligible. This may introduce 
significant simplifications in near-wall modeling, where 
this limit applies as y + — * 0. 



F. Relation to Rotta [4] 

The present result in (56) with (57) for ReA > 1 or 
(58) for ReA — ^ is the first nonlocal pressure-strain 
correlation that rigorously accounts for the effect of spa- 
tial variations in the mean velocity gradients on the tur- 
bulence anisotropy. Previously, Rotta [4] derived some 
of the components of M^ k and mp M^ fc , though not 
enough to construct even the leading nonlocal term in 
(50). In particular, Rotta used an inertial-range form 
for f(r) similar to (35) and the Gaussian form in (41) 



E 

Tl=2 



r W-q 

Ly 2 O i j 



(59) 



C 



{n) 



C x - 



In homogeneous flows, for which prior purely local mod- 
els for Uij have been relatively successful, the Lapla- 
cians of Sij in (59) vanish, and thus the present non- 
local pressure-strain formulation recovers the local form 
in (21), since C| X) = 4/5 in both (57) and (58). For inho- 
mogeneous flows, when (59) is introduced in (2) it gives 
a new anisotropy transport equation that accounts for 
both local and nonlocal effects via the present fundamen- 
tal treatment of spatial variations in the mean velocity 
gradients in (10). Note in (2) that the definition of Pij 
withP = P rm /2: 



p. 



u^u'j gives 
= -^S~j + (aaWij - W u aij) (60) 



auSij + Suaij 
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where the mean- flow rotation rate tensor Wij is given by 



— 1 / dui du, 

w,. - ' 



2 \ dxj dxi 



(61) 



From (60), the production terms in (2) thus require no 
additional closure modeling, while current standard mod- 
els summarized in Refs. [1, 2] may be used for the remain- 
ing eij and Dij terms. 

However, (59) does not account for possible additional 
anisotropic effects in since the present nonlocal 

pressure-strain result in (56) is based on the central hy- 
pothesis that Ru(r) in (26) can be represented by its 
isotropic form in (18). Fundamentally-based approaches 
for any such remaining anisotropic effects in (59) have yet 
to be rigorously formulated, however it has been heuristi- 
cally argued (e.g., [9, 19]) that such additional anisotropy 
effects may be represented by higher-order tensorial com- 
binations of ciij, Sij, and W^. The most general of such 
combinations that remains linear in o^- is 



aniso) 



k » 



= < ' ; I auSij + Siiaij - -a n iS n iSij^j (62) 
+C 4 (auWij - Wuaij) , 



where the constants C 3 and C4 can be chosen to pre- 
sumably account for such additional anisotropy effects. 
In general, choices for these coefficients vary widely from 
one model to another; a summary of various such models 
is given in Ref. [2]. 

When (59) is combined with (62), it provides an 
anisotropy transport equation that accounts for both lo- 
cal and nonlocal effects, as well as possible additional 
anisotropy effects, in the pressure-strain correlation as 



Daij e - 

—- = ^ ai - k a t]+ a 2 S. 



(63) 



E 

n=2 



a 



(n) 



(v 2 )' 



2 " 






+ a 3 ^ 



-a 4 (auWij - Wuaij) + - 



2 

Dij — ( ciij + —Sij 



D 



where the coefficients are given in (57) or (58), and 
the cti arc defined as 



a\ = — - 1 + C\ , a 2 = C 2 

otz = C3 - 1 , a 4 = C 4 - 



(i) 



4 

3 ' 



(64) 



1. 



Values for the constants C\, C3 and C4 in (64) may be 
inferred from prior purely local models, such as the Laun- 
der, Reece and Rodi (LRR) [19] or Speziale, Sarkar and 
Gatski (SSG) [9] models, which are all based on forms of 
(63) without the nonlocal effects given by the series term. 
However, optimal values for these constants may change 



in the presence of the nonlocal pressure-strain term in 
(64). 

With respect to the remaining terms in (63), for high 
Reynolds numbers the dissipation tensor is concen- 
trated at the smallest scales of the flow, which are as- 
sumed to be isotropic. Thus, consistent with the central 
hypothesis on which the present result for the pressure- 
strain tensor is derived, the dissipation is commonly rep- 
resented by its isotropic form = ( e -5-> [2, 6]), 
with the result that the dissipation term in (63) vanishes 
entirely. The only remaining unclosed terms when (63) 
is used with the ensemble-averaged Navier-Stokes equa- 
tions are the transport terms and D, and these are 
typically represented using gradient-transport hypothe- 
ses, with several possible such formulations summarized 
in Ref. [2]. 

A number of different approaches can be taken for solv- 
ing (63). First, this may be solved as a set of six coupled 
partial differential equations, together with the ensemble- 
averaged Navier-Stokes equations, to obtain a new non- 
local Reynolds stress transport closure that improves on 
existing approaches such as the LRR and SSG models 
in strongly inhomogeneous flows. Alternatively, equilib- 
rium approximations may be used to neglect the Daij/Dt 
and terms in (63) to obtain a new explicit nonlocal 
equilibrium stress model for a^, analogous to the exist- 
ing local models developed, for example, by Gatski and 
Speziale [20], Girimaji [21], and Wallin and Johannson 
[22] . Perhaps preferably, a new explicit nonlocal noncqui- 
librium stress model for can be obtained from (63) 
following the approach in Ref. [23], by explicitly solving 
the quasi- linear form of (63), namely 



Da^ e - 

~5T = ~ ai k a v +a ^ S v 



(65) 



+ E 

n=2 



a 



(n) 



(v 2 )' 



In so doing it is possible to obtain a new explicit form 
for the anisotropy that accounts for both nonlocal 
and nonequilibrium effects in turbulent flows. 



IV. CONCLUSIONS 

A new rigorous and complete formulation for the rapid 
pressure-strain correlation, including both local and non- 
local effects, has been obtained in (56) with (57) for 
Re\ > 1 or (58) for Re A — > 0, and with A in (37) 
and (40). Nonlocal effects are rigorously accounted for 
through Taylor expansion of the mean velocity gradients 
appearing in the exact integral relation for in (10). 
The derivation is based on the central hypothesis that the 

(r) 

nonlocality in 11^ ■ is substantially due to spatial varia- 
tions in duk/dxi in (10), and that in order to address 
this effect all other factors in (10) can be adequately rep- 
resented by their homogeneous isotropic forms. The re- 
sulting rapid pressure-strain correlation in (56) takes the 
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form of an infinite series of increasing-order Laplacians of 
the mean strain rate field Sij (x), with the n = 1 term re- 
covering the classical purely-local form in (21), and with 
the remaining n > 2 terms accounting for all nonlocal 
effects due to spatial variations in the mean-flow velocity 
gradients duk/dxi. 

Aside from the central hypothesis on which the present 
approach is based, the sole approximation lies in the need 
to specify a form for the longitudinal correlation function 
fir). The particular specification does not affect the fun- 
damental result in (56), and serves only to determine the 
pressure-strain coefficients For the classical expo- 

nential form in (39) appropriate for Re a ^> 1, the corre- 
sponding coefficients are given in (57), while for the exact 
Gaussian form in (41) appropriate for Re\ — > the co- 
efficients are given in (58). The integral scale A in (56) 
determines the size of the region around any point over 

(r) 

which nonlocal effects are significant in IIL . In general, 
A can be obtained via (37), with C\ in (40) giving good 
agreement with the inertial-range form of fir) in (35) 
and (36) for Re A > 1. 

The agreement of the present n = 1 term with the 
purely local form in (21) obtained by Crow [5], and 
with the limited components obtained by Rotta [4] for 
the leading (n = 2) nonlocal term, support the valid- 
ity of the present derivation. The present results, how- 
ever, go much further by accounting for all components 

(mp...)^iijk f° r au n i which together have allowed the 
complete form of both the local and nonlocal parts of 



the rapid pressure-strain correlation ILL to be obtained, 
within the central hypothesis on which the present ap- 
proach is based. The present result thus gives the first 
rigorous nonlocal form of the rapid pressure-strain corre- 

(r) 

lation IIL for spatially varying mean velocity gradients 
in turbulent flows. 

Using the present result for n^' in (56) with (37) and 
(40) and with (57) or (58), a nonlocal transport equa- 
tion for the turbulence anisotropy has been obtained in 
(63) and (64). The resulting nonlocal anisotropy equa- 
tion can be solved by any number of standard meth- 
ods, including full Reynolds stress transport closure ap- 
proaches, algebraic stress approaches, or the nonequi- 
librium anisotropy approach outlined in Ref. [23] based 
on (65). This nonlocal anisotropy equation should give 
significantly greater accuracy in simulations of inhomo- 
geneous turbulent flows, including free shear and wall- 
bounded flows, where strongly nonuniform mean flow 
properties and significant large scale structures will in- 
troduce substantial nonlocal effects in the turbulence 
anisotropy. 
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